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The combination of density functional theory with other approaches to the many-electron prob- 
lem through the separation of the electron-electron interaction into a short-range and a long-range 
contribution is a promising method, which is raising more and more interest in recent years. In this 
work some properties of the corresponding correlation energy functionals are derived by studying 
the electron-electron coalescence condition for a modified (long-range-only) interaction. A general 
relation for the on-top (zero electron-electron distance) pair density is derived, and its usefulness is 
discussed with some examples. For the special case of the uniform electron gas, a simple param- 
eterization of the on-top pair density for a long-range only interaction is presented and supported 
by calculations within the "extended Overhauser model". The results of this work can be used to 
build self-interaction corrected short-range correlation energy functionals. 
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I. INTRODUCTION 



In recent years, there has been a growing interest in ap- 
proaches that combine density functional theory Q, |j| 
(DFT) with other approximate methods to treat the 
many-electron problem. In most cases, this combination 
is achieved by splitting the Coulomb electron-electron in- 
teraction 1 / ri2 into a short-rEinge (SR)and a long-range 
(LR) part (see, e.g., Refs S H H i, Ho, H, ll- 
The idea is to use different, appropriate approximations 
for the long-range and the short-range contributions to 
the exchange and/or correlation energy density function- 
als of the Kohn-Sham (KS) scheme, to treat, e.g., near- 
degeneracy effects or van der Waals forces. These ap- 
proaches are often inspired by the observation that long- 
range correlations are not well treated by local or semilo- 
cal density functionals, but can be dealt with by other 
techniques, like the standard wavefunction methods of 
quantum chemistry. Conversely, correlation effects due 
to the short-range part of the electron-electron interac- 
tion can be well described by local or semilocal function- 
als (appropriately modified). 

The error function and its complement (see Fig. 
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are often used 

H II i, in im for the splitting of the 
Coulomb interaction, since they yield analytic matrix el- 
ements for both Gaussians and plane waves, i.e., the most 
common basis functions in quantum chemistry and solid 
state physics, respectively. The parameter fi controls the 
range of the decomposition. Correspondingly, the uni- 
versal Coulombic functional of the electron density n(r) , 
F[n], as defined in the constrained search formalism |l3| . 



F[n] = min(^'|T + V;e|*) 
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FIG. 1: The splitting of the Coulomb interaction l/ri2 into 
a short-range (SR) and a long-range (LR) part as defined in 
Eq. CJ. Here /i = 1. When ^ ^ oo we have «sr — > and 
ULR l/ri2, and when fj, —> we have vsr —> l/ri2 and 
"LR 0. 



tary short-range part, F[n] = F^p [n] + Fg^^[n], 
Fi:^[n] = min (vI/^'|T + ^^1,1*^) 



FskM = F[n]-Fi:^[n], 



(3) 



or, alternatively, into a short-range part and a comple- 
mentary long-range part 



mm 
F[n] 
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can be divided into a long-range part and a complemen- 



These two decompositions lead to different exchange- 
correlation energy functionals that need to be approxi- 
mated; they are compared in Ref. '14| where their advan- 
tages and disadvantages are discussed. 

In the present work we focus on the properties of the 
long-range and short-range correlation functionals that 
come from the modification of the electron-electron in- 
teraction at short distances, i.e. those properties that 
are due to the change in the electron-electron coalescence 
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conditions. This means that we are only concerned with 
the functional of the decomposition of Eq. (|3J), which 
involve a many-body wavefunction ^E*^ of a system with 
an electron-electron interaction that is softer than l/ri2 
for small ri2 (see Fig.0. This decomposition is the one 
used in the approaches of Refs. IM rTllsllfll I1(TI IllL 

The paper is organized as follows. In Sec.|n]we define 
the long-range and short-range correlation energy func- 
tional that are the object of the present study. In Sec. lIIII 
we analyze the short-range properties of the pair den- 
sity of a general many-electron system with interaction 
erf(/xri2)/ri2 when gets larger and larger: we derive 
an expansion for ^ — > oo of the on-top (zero electron- 
electron distance) pair density, and, following (and partly 
correcting) the work of Ref. lu, an expansion for /i — > oo 
of the correlation energy functionals. Section IIVI is de- 
voted to the special case of the uniform electron gas: 
starting from the exact high-density limit, a simple pa- 
rameterization of the on-top pair density as a function 
of fi is proposed, and is favorably compared with the 
results obtained from the "extended Overhauser model" 
^^01 for the same quantity. The last Sec. Ivl explains, 
with some examples, how the results of this work can be 
used to build self-interaction corrected approximations 
for short-range correlation functionals. 

Hartree atomic units are used throughout this work. 



II. DEFINITIONS AND BASIC EQUATIONS 

When the universal functional F[n\ is decomposed as 
in Eq. (PJ, we have a model system, whose wavefunc- 
tion is denoted 5"^, which has the same density n(r) 
of the physical system and electron-electron interaction 
erf(^ri2)/ri2. When — > this model system becomes 
the Kohn-Sham system, with no electron-electron inter- 
action, while when /i cx) the model system approaches 
the physical one, with interaction l/ri2. By definition, 
the density is the same for all values of /u. In the ap- 
proach of Refs. Qj IS Hi ^ the model system at a fixed 
fi is treated with a multideterminantal wavefunction. In 
general, if /i is not too large, few determinants describe 
^''^ quite accurately (because of the smaller interaction, 
and also because of the absence of the electron-electron 
cusp); the larger is the chosen value of /i, the larger is the 
needed configuration space and thus the computational 
cost. The remaining part of the energy is provided by 
the complementary functional FgR[n] = — F-[^^[n] 
of Eq. (PJ, which can be divided into Hartree, exchange, 
and correlation contributions in the usual way (for an 
alternative separation of exchange and correlation, see 
Ref. 113), 

^h,srW - IJ'^'^J dr'n{r)nir')vti^i\r-r'\) (5) 

KsrW = ('i>l^sRl'i'>-^H,sRW (6) 

-^c,srW = ■P'srW - -E^h.srW - -E^L'^srW' (7) 



where $ is the Kohn-Sham determinant. Notice that the 
Hartree and the exchange functional are linear in the in- 
teraction, so that -Eh.srN — srN ^^'^ ^a-, Sr["] = 
^^,SrW' where ^^^^^srN and -B^^srW are the short- 
range functionals of the decomposition of Eq. Q). The 
correlation energy, instead, depends on the wavefunc- 
tion ^''^ and we thus have e'^ g^ [n] ^ -^'cSrI'^]- The 
complementary correlation functional E'^ srI"'] the dif- 
ference between the usual Coulombic correlation energy 
-Bci?^], and the long-range correlation energy functional 

KlrM = mr +¥^^1^^^) - mr +v^^m (s) 

KsrW = EM^Kli^H (9) 

In what follows we study the short-range functional 
^cSrM or, equivalently the long-range functional 
-^c lrI*^]' the limit of large /x, i.e., when the model 
system described by is approaching the physical sys- 
tem. 

Following Toulouse et al. we start from the 

Helmann-Feynmann theorem which gives 

Cm V^r Jo 

(10) 

where the spherically and system-averaged pair density 
(or intracule density) f^{r 12) is obtained by integrating 
l^fp over all variables but ri2 — |r2 — ri |: we first define 
the spherical average of the pair density (with r = ri), 

j |*^(r,ri2,r3,...,rAr)|2^dr3...dr^, (11) 

and then integrate over all reference positions r, 

r(ri2) = j P^{v,r,2)dv. (12) 

The correlated part of f^{r 12) appearing in Eq. lfTU|l is 
fc = f - /ks, where /ks(?'12) is obtained by replacing 
^E''^ with the Kohn-Sham determinant $ in Eq. pi|l . 

The correlated intracule f^{r 12) can be expanded in 
its Taylor series around ri2 = up to some order M, 

M-l 

/en'^12)= E ^»(^) ^12 +0(^*2)- (13) 

By inserting this expansion into Eq. 11U|) we find 




(14) 

This means that when /i — > cx3, i.e., when we are ap- 
proaching the full interaction l/ri2, the correlation en- 
ergy functional -E^ srM ^^s an expansion in powers of 
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/Lt~^ whose coefficients are determined by the short-range 
behavior of In order to determine as many coeffi- 
cients as possible in the expansion of Eq. p4|l we thus 
have to know how the c„ (/i) behave in the hmit of large 
11. In Ref. ^3 the same expansion was considered, and 
the first two terms were obtained by simply inserting in 
Eq. (|14|l the values of cq and ci for the physical system 
(with Coulomb interaction). In the next Sec. lIIII we study 
the general problem of determining the short-range be- 
havior of f^iry^) in the limit /Lt — *■ oo. We find the same 
result of Ref. for the first term (cx ^^^) of Eq. 
but a different result for the second term (cx fJ.~^), and we 
explain why. Moreover, we obtain the first-order (in 
term of the expansion for large /i of the on-top /''(O). 



III. SHORT-RANGE BEHAVIOR OF A SYSTEM 
WITH INTERACTION eri {/ir 12) /r 12 WHEN ^ ^ 00 

We start from the Schrodinger equation for the wave- 
function 5*^, 



Schrodinger-like equation reads 



(P erf(/Lia;) 



u^(x) =£^(x,R,r3,...,r7v)ii'^(a;), 

(16) 

where £^ is a complicated operator that does not affect 
the result as long as it remains bounded when /i — > 00 and 
X — > 0, as it is reasonable to assume from the hamiltonian 
of Eq. 1113 We change variable y = fix, and 

divide both members of Eq. H16|l by fi^ to obtain 
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u''(2/) = — £^(y,R,r3,...,rjv)u''(y). 

(17) 

We then expand u'^''{y) for large /i, 

u^iy) - u^°^\y) + j^u^-'Hy) + O (^) , (18) 

insert this expansion into Eq. H17|l . and impose that the 
left-hand-side be of order fi~^ as the right-hand side. 
With the boundary condition that ip'^{x) is finite at 
x = 0, we obtain 



(15) 
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where the one-body potential ^^^(r) keeps the density 
equal to the one of the physical system for every fi. 

When fi ^ 00 the interaction erf(/xri2)/ri2 gets larger 
and larger for small ri2 (ri2 <C l//x) . If we thus fix a finite 
but very large value of /x, we can use the same arguments 
that lead to the derivation of the electron-electron cusp 
condition for the Coulomb interaction [3 0, 120I l2lj| . 
i.e., we can isolate two coalescing electrons (say, 1 and 
2) in the hamiltonian of Eq. H15|l . and switch to vari- 
ables ri2 = ri — r2 and R = (ri -I- r2)/2. In the limit 
^12 = |ri2| — > there must be a term in H^^'^^ that 
compensates the divergence (or, more precisely, the very 
large value) of erf(/iri2)/ri2. As for the Coulombic sys- 
tems, this compensation comes from the relative kinetic 
energy term, and to determine the small ri2 behavior of 
the spherical average of jvp'^p we only need to look at 
the Schrodin ger equation for the relative motion of two 
electrons [l^ [l^ [2(J| approaching each other with rela- 
tive angular momentum i — Q. Higher ^, in fact, will 
contribute to |\E'^p to orders in the limit of small ri2. 
Only in the case of a fully polarized system the case 1=1 
must be considered to determine the ri2 — *■ behavior 
of l^'^^p, since only odd I are allowed |0| (triplet sym- 
metry); this case is considered in Appendix 151 The rare 
case of unnatural parity singlet states [23 (which needs 
£ = 2) is not considered in this work. 

As it was done for the Coulomb electron-electron in- 
teraction 0, 0, [23 , we thus focus on the relative wave- 
function ip^{ri2) for two electrons in the t = state. 
By defining x = ri2 and u^{x) = xij}^{x), the relevant 
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d^u^-^^y) 
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aerf(y). 



(19) 
(20) 



and we find that the final solution for ^'^{x) from 
Eqs. (113-1201) is 



i;^'{x) 



1 + xpi{fix) 



7^ + - 



(21) 



where Ai is a constant coming from the integration of 
Eq. (|20|l that is not determined by the condition that 
ip^{x) is finite in a; = 0, and a determines the value V'(O) 
for the Coulombic system (fi = 00). The function pi{y) 
is given by 



Pi (2/) 



e-y - 2 



1 



2 ; erf(2/). 



2^y ■ V2 ■ 4y 
and has the following asymptotic behaviors 



pi(y^o) = 
pi{y ^ 00) = 
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(23) 
(24) 



The spherically- and system-averaged pair-density of 
Eq. H12|) has thus, to leading orders in for large /z, 
the small-r'12 expansion 



r(ri2) = /(O) 



1 -I- 2ri2Pi(^ri2) 



2 2Ai 



(25) 

where /(O) is the on-top value corresponding to the full 
interacting system [/(O) is proportional to a^, where a 
determines u^°°\y) in Eq. lfTO|l ]. Equation (|^ tells us 
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that if in Eq. ^ we fix ri2 equal to a small value tq ^ 1, 
and then let /i go to oo we recover the Coulombic cusp, 
/(''o) = /(0)(1 + ^0 + •••)■ But for any finite large /z, 
Eq. (|23|l shows that we always obtain a quadratic behav- 
ior for small na, /^(ria) = /(0)(1 + 2rf^^i/3^ + ...). 
This is how the cuspless wavefunction corresponding 
to the interaction evi {fir 12) / r 12 develops the Coulombic 
cusp in the — > 00 limit. An alternative derivation 
of Eq. (|25() . more similar to what one usually does for 
the Coulombic cusp HI El US EH, is reported in Ap- 
pendix^ 

To obtain the complementary short-range correlation 
functional we can insert Eq. H25I) into Eq. (I14II . which 
gives 
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(26) 

We see that to fully determine the term cx /i ''in Eq. H26() 
we have to know the value of the constant Ai. This 
constant determines how the on-top /^(O) approaches the 
Coulombic value /(O) for large jjL. In fact, from Eq. 
we have 



(27) 



In Ref.T^ the fji — > 00 limit of the long-range interaction 
was formally rewritten as the Coulomb interaction l/ri2 
plus a perturbation [Tol | 



= 2S{ri2) + O — 

ri2 ri2 fi^ Vm 



(28) 



where 5{r) is the Dirac delta function. The fact that the 
perturbation is of order lead to the conclusion 
that the perturbation on VE"^ is also of order with 
respect to the Coulombic case, which would correspond 
to Ai = —l/y/ir in Eq. I|27|l . However, because of the 
singular nature of the Dirac delta function, this argument 
does not hold at ri2 = 0. 

To determine the correction to the on-top value when 
fi 00, here we take a large value of fi and a small value 
ri2 — ro such that //^' ^ ro ^ 1 (take, e.g. tq — l//i'~'^ 
with < g < 1). For such value of tq we have /iro 3> 1 
so that from Eqs. ()24|l and we obtain 



riro)=f{0) 1 -f ro + 



2Ai 



(29) 



where the next leading terms are of order and/or 
Tq. We then notice that Ai cannot be equal to — l/y^, 
since any value of < would make f^{ra) smaller 
than the full interacting value f{ro), while, because for 
small ri2 the long-range interaction erf(/iri2)/ri2 is less 
repulsive than l/ri2, for ro small enough we expect that 
f^ifo) > /(''o)- So Ai > 0, and thus the correction to 
the on-top value must be of order However, the 

argument of Ref. should be valid when ri2 » 1//^. 
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FIG. 2: The Coulombic cusp /(0)(1 + r^i) (with /(O) = 1) 
is compared to the expansion of j'^irrz) in the /i 00 limit, 
F{rxi) = /(0)(1 2ri2Pi(/iri2) + 2/(0F^i) + 2^i//x), with 
A\ = 0. Any value Ai > makes j^(r\2) differ from the 
Coulombic jirx-i) of orders also in the region ri2 ^ 



That is, we still expect from Eq. I|28|l that the perturbed 
differs from the Coulombic 'F of an order higher than 
for rx2 ^ l/^i- This is achieved only if Ax — 0, as 
shown by Eq. H29() . The result corresponding to j4i = 
is illustrated in Fig. |21 where we compare the Coulomb 
cusp /(0)(l + ri2) to the short-range expansion of /^(ri2) 
of Eq. (|25|1 . with A\ — 0. Any value of A\ larger than 
zero makes the difference between the Coulombic f{r 12) 
and /^(ri2) of order also in the region ri2 ^ 
while with Ai = (as in Fig. [21) this difference is of higher 
order, as expected from Eq. H28|l. 
We thus conclude that 



no) = /(o) 1 + 



'TT fl 



(30) 



This equation is also confirmed in the next Sec. II VI for the 
case of the high-density electron gas that can be treated 
exactly. 

The final expansion of the short-range functional 
sr['^] fo^' large fj, is then 
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n] = /c(0)^ + /(0) 
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(31) 



where /(O) and /c(0) are the on-top value and its corre- 
lated part, fc — f — /kSj of the physical system. This 
expansion differs from the one of Ref. 10 by a factor \/2 
in the second term (see Appendix ^ for comments on 
this discrepancy). The two expansions for the case of 
the He atom are compared in Fig. 13 with the "exact" 
results [i[l3 forl'^ sRN = ^^--E^.lrN- The "exact" 
on-top value /(O) is taken from Ref. 22. We see that 
the new expansion more accurately reproduce the "ex- 
act" data for large n, and that the /j,-value for which it 
breaks down (i.e., where the expansion has its minimum) 
coincides with the minimum value of /i for which it still 
gives accurate short-range correlation energies. 
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FIG. 3: The "exact" values of -Bc.srN = Ec - £^c,lrW 
the He atom pLlinl| are compared with the large- /i expansion 
of Eq. II31II . and with the previous result for the same ex- 
pansion given in Eq. (30) of Ref. ITol In both expansions the 
"exact" on-top value is taken from Ref. 



In general the on-top value of the physical system, 
/(O), is not accessible. A plausible approximation pro- 
posed in Ref. |2 consists in replacing /(O) in Eq. 
with its local-density approximation (LDA) value, 



/lda(O) 



n{vf g{ri2 = 0;7i(r))dr, 



(32) 



where g{r i2\n) is the pair distribution function Il6l l23|| 
of the standard electron gas model (with Coulomb in- 
teraction l/ri2) of uniform density n. The new Eq. H30(l 
allows to estimate the physical on-top value starting from 
the one of the model system Potential applications 
of this idea are discussed in Sec.0 together with simple 
examples. Notice also that Eq. H3U|) is also valid locally, 
i.e., we have 



P2^(r,ri2 =0) = F2(r,ri2 =0) 1 



/tt /i 



(33) 

where Pl^iv^r 12) was defined in Eq. Ullfl . and P2(r,ri2) 
is the pair density of the physical system (/i = 00). 

In Appendix El we also consider the case of a fully 
polarized system, for which we find that the short-range 
correlation energy has the large-/z expansion 
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where /"(O) and /"(O) are the second derivative at ri2 = 
and its correlated part of the physical, fully interacting, 
system We also found that, again only in the case of 
a fully polarized system, the second derivative of f^{r 12) 
at ri2 = approaches the one of the Coulombic system 
as 



(/T(o)^/"(o)(i + ^) + 
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(35) 



IV. ON-TOP PAIR DENSITY OF A UNIFORM 
SYSTEM WITH INTERACTION erf {fir 12) /r 12 

Before coming to applications, we consider the special 
case of the uniform electron gas, for which something 
more than Eq. (|30|l can be done. We consider a uniform 
system with long-range-only interaction. 
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(36) 



where is the modified electron-electron interaction 



(37) 



is, accordingly, the interaction between the electrons 
and a rigid, positive, uniform background of density n = 
(4^r3/3)-i 



N 



erf(^|ri - x|) 



(38) 



and Vj^ is the corresponding background-background in- 
teraction 

When /J, — > 00 tends to the standard jcUium hamil- 
tonian, while when /i — > we recover the noninteracting 
electron gas. 

We focus on the /i-dependence of the on-top value of 
the pair-distribution function |23| g{ri2 = 0, rs,/x), 
which has its own interest to construct the LDA approx- 
imation for the long-range and short-range functionals, 
and for spin-density functional theory in the framework 
of the alternative on-top interpretation |25| . The relation 
between the function g and the function / of Eq. I|12() is 
g = 2f/nN. 



A. High-density limit 

As in the Coulomb gas, by switching to scaled units 
Si = Vi/vs, we see that when — > the potential of 
Eqs. (|37|) - (|39|l becomes a perturbation to the noninter- 
acting gas. Defining the correlation contribution to the 
on-top value, gc{0,rs,n) = (?(0, rs,/Lt) — i, and following 
Kimball j23|, the first-order correction (with respect to 
the interaction potential) to the on-top pair density is 



5c(0, 
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0,/^) = 6 / ASDit,fi)t^ 
Jo 



(40) 



where ASd (t, ^) is the direct second-order contribution 
to the static structure factor 26], 
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FIG. 4: The function h{z) that determines the high-density 
limit of the on-top pair density of the "erf" gas [see Eq. I|42|l ] . 
The numerical evaluation of Eq. 14011 (points) is compared to 
the fitting function of Eq. 1451 (solid line) . 



np is the usual Fermi occupation function |2a | and Wee (q) 
is the Fourier transform of the electron-electron inter- 
action. In Eq. H4Q|I the scaled variable t — q/2kF 
[kp = {ars)~^, a = (^)^/"^] has been used. The function 
ASdUjU) is thus equal to the one computed by Kim- 
ball [2g and reported in his Eq. (11), except for a mul- 

+ 272/2. 

tiplying factor e~ "^f/^ coming from the Fourier trans- 
form of erf(/iri2)/ri2. From Eqs. (|40|l and (|41() we find 



(42) 



where the function h{z) has the following asymptotic be- 
haviors 



h{z ^ 0) 
h{z oo) 



-^(1 -ln2)z 
a 
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0{z') (43) 
+ 0{z-^), (44) 



and flHD = -a(7r2 61n2 - 3)/57r w -0.36583 is the 
high-density limit of the standard jellium model. Notice 
that Eq. (|44|l confirms, for the high-density electron gas, 
Eq. 1201). 

For intermediate values of z we computed numerically 
the function h{z)^ and we found that it can be accurately 
fitted by the Fade form 



h{z) 



ai z^ + 02 



l + biz + b2z'^ + b3Z^' 



(45) 



with ai = -^(1 - In 2), 6i = 1.4919, feg = 1.91528, 
0-2 = ciHD^3, ^2 = (fli ^ b^cv / ^/tt) / a-HD . The numerical 
results and the fitting function of Eq. (|45|l are reported 
in Fig. 21 



B. Interpolation formula 

The high-density hmit of Eq. and of Fig.^tells us 

how (at least for small r^) f/c(0, Ts, fi) approaches the two 



FIG. 5: The correlation on-top pair density of the electron 
gas with interaction erf (/xri2)/ri2 for different /i (in a.u.) as 
a function of the dimensionless density parameter r^. The 
results from the Overhauser model are compared with the 
interpolation formula of Eq. 1)47^ . The dotted line corresponds 
to the standard jellium model with interaction l/ri2. 



limits, the noninteracting gas (/i — > 0) and the Coulomb 
gas (/i oo). 

A simple interpolation formula for all densities can be 
built by assuming that the dependence of gc{0,rg, fi) 
is roughly the same at each r^. We thus start from the 
parameterization of the on-top pair density of the jellium 
model given in Rcf. 13 



1 



g{0, rs, ^ = oo) -(1 - Br^ 



i?r^)e-*% 
(46) 

where C = 0.08193, D = -0.01277, E = 0.001859, 
d = 0.7524 and B ~ —2aHD — d, and we simply 
rescale homogeneously all the coefficients with the func- 
tion h{z — /i//ci?)/aHD, 



=(0,rs,M) 



-drs /i(z)/aHD p 



1 



Hz) 



"HD 



-D 



h{zf 



''HD 



'■HD 



(47) 



This simple guess smoothly interpolates between the fx 
and — > oo limits, and is exact when — > 0. 



C. Results from the Overhauser model 

To check the validity of the interpolation formula of 
Eq. (|47|l we evaluated the on-top Of,(0 , r.^. li) within the 
"extended Overhauser model" [ia)l3> which gave good 
results for the standard jellium model and for two- 
electron atoms Notice that the on-top value is not 
known exactly. The differences between the jellium on- 
top pair densities from different approximate methods 
(including the "extended Overhauser model") are dis- 
cussed in Refs. m H IS 

The scattering equations of the "extended Overhauser 
model" are widely explained in Refs. ll^lsil Here we sim- 



7 



ply solved the same equations with the electron-electron 
interaction evf { fj,r 12) /r 12 screened by a sphere of radius 
of uniform positive charge density n and attracting the 
electrons with the same modified interaction, 



Veff(ri2,rs,/z) 



erf (/xri 
ri2 



^erfMr--ri,|)^^,^ 
r' - ri2 



This potential is reported in the Appendix of Ref. l27l 
where it has been used for two-electron atoms with rather 
accurate results for the corresponding short-range corre- 
lation energy. Kff (ri2, Ts, /i) is a screened potential that 
tends to the "Overhauser potential" when ji ^ 00, 

and which goes to zero when /i ^ 0. As in the original 
Overhauser model, the idea behind Eq. H48|l is that the 
radius of the screening "hole" is exactly equal to . 

The results for the on-top (7c(0, rs,/x) from the Over- 
hauser model are reported in Fig. as a function of 
for different values of /i. We see that the the simple in- 
terpolation formula of Eq. (|47|l accurately captures the 
/i and Vs dependence of gc(0, Tg, 



APPLICATIONS, PERSPECTIVES, AND 
CONCLUSIONS 



The main results of this work are (i) the corrected ex- 
pansion of the short-range correlation energy functional, 
Eq. H31|l . (ii) the expansion of the on-top pair density 
of Eqs. (|30|l and H33|) . and (iii) the parameterization of 
the /Lt- dependence of the on-top pair density of the uni- 
form electron gas, Eq. H47|l . All these results (i-iii) can 
be useful for the construction of approximate short-range 
correlation e nerg y functionals: 

(i) In Ref. [13 the large-/i expansion of the correlation 
energy functional has been already used to construct ap- 
proximations: the idea is [Tol | to interpolate between a 
given density functional approximation (DFA) of stan- 
dard KS theory [s^l at /i = 0, and the /i — > 00 expansion 
of -Ef, srM- III the spirit of the usual DFT approxima- 
tions, this interpolation is done locally, i.e.. 



c,sm 



i] = y drn{r) 



SR 



(r), 



where sr('") built, e.g., as [lol | 



,DFA 



^c, SR 



1 + di fi + d2 fi'^ 



(49) 



(50) 



The parameters di and ^2 are fixed by the condition 
that Eq. (|49|) recovers the correct /i — > 00 expansion 
of srWi ^^'^ eJJ'^^ can be, e.g., the PBE correla- 
tion functional js^ of standard Kohn-Sham theory or 
any other available approximation. This correlation func- 
tional can be combined with a similar interpolation for 
exchange ^3 : o'' with the exchange functional of Heyd, 
Scuseria, and Ernzerhof 6] . This way of constructing ap- 
proximations should be improved by using the corrected 
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FIG. 6; The "exact" values of -Bc.srN = Ec ~ -B^_lrN for 
the Be atom pllT(i^ are compared with the large-p expansion 
of Eq. 13111 . and with the previous result for the same ex- 
pansion given in Eq. (30) of Ref. ITol In both expansions the 
"exact" on-top vahie is taken from Ref. 



expansion of Eq. , as suggested by Fig. |3| In Fig. 
we also show similar data for the Be atom: again, the 
corrected expansion is closer to the "exact" data 1^ IhI 
at large fi than the previous expansion used in Ref. Uoi 

(ii) To impose the correct large- /Lt expansion in approx- 
imations like the one of Eq. H5()|l we need an estimate of 
the physical (/x = 00) on-top pair density P2(r,ri2 = 0). 
In Ref. 10 the LDA approximation [the integrand of the 
right-hand side of Eq. H32() ] was used. The new Eq. H33() 
allows to use the on-top pair density of the partially 
correlated wavefunction '^'^ to estimate P2 (r, ri2 = 0) . 
In fact, once we have made a calculation at a given 
(moderately large) /i (say, /i — /xq) we can estimate 
^2(r,ri2 = 0) as 



P2(r,ri2 =0) 



Fr(r,ri2 = 0) 



(51) 



There are cases in which this estimate could be much bet- 
ter than the one obtained by using the LDA approxima- 
tion for the physical on-top pair density. In fact, the use 
of the partially correlated P^(r,ri2 = 0) would correct 
the self-interaction error of LDA, becoming exactly equal 
to zero for any one-electron density. Consider the exam- 
ple of the stretched II2 molecule, for which the estimate 
from Eq. (|51|l would be essentially exact (equal to zero 
for any /i > 0), while LDA gives a spurious nonzero on- 
top value, unless we consider the spin broken-symmetry 
solution. 

To show that Eq. (|51|l gives indeed a quantitative re- 
liable estimate of the physical on-top pair density we 
have considered the simple example of the He atom, and 
we have inserted Eq. (|51l) directly in the expansion of 
Eq. H31|l : the error on the estimated short-range correla- 
tion energy at /.to = 2.5 is 3 mH; if we choose /.tg = 2 the 
error in ^^c^brN ^ txiR, and for /xq = 1 is H rnH. Of 
course, when /iq becomes too small, the large-/z expansion 
of Eqs. (|31|l and H51|) is no longer valid. 
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(iii) The on-top pair density of the uniform electron 
gas with long-range- only interaction of Eq. (|47|l can be 
used, in combination with the correlation energy of the 
spin-polarized long-range electron gas '3 4*1 , to implement 
the local approximation for the on-top pair density inter- 
pretation of spin density functional theory [25j |. 

Another interesting application of Eq. (|47|l is con- 
nected to point (ii): the construction of functionals that 
explicitly depend on the on-top /^(O) [or locally on 
P2'(r, fi2 — 0)] of the partially correlated wavefunction 
could use the p. dependence of the on-top LDA value to 
go beyond Eq. (|^ . 



/c(0) 



^ g^^Q' ""'^^^^ ii = oa)dv 



For example for the He atom Eq. 1)521) at /io = 2 gives 
/^(O) = -0.086, while Eq^jnUl gives -0.090. The corre- 
sponding "exact" value [22] is —0.085. Local versions of 
Eq. ifS^ can be also considered, e.g.. 



/,(0)« / n(r)^Pf,(r,ri2=0) 



5c(0,rs(r),/^o) 



dr, 

(53) 

where P^°(r, ri2) is obtained by subtracting the Kohn- 

Sham pair density from . Again, the advantage of 
including in the construction of short-range functionals 
the on-top P^^(r,ri2 = 0) is to locally remove the self- 
interaction error. 

In conclusions, we have presented a comprehensive 
study of the short-range behavior of systems interact- 
ing with the potential erf(/iri2)/ri2, in connection with 
the properties of long- and short-range correlation en- 
ergy density functionals. The same kind of analysis can 
be of course repeated for other splittings of the Coulomb 
electron-electron interaction 0, l35j| . Future work will be 
mainly oriented to the exploration of the promising ap- 
proach of short-range functionals that explictly depend 
on P2^(r,ri2 = 0). 
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APPENDIX A: ALTERNATIVE DERIVATION OF 
EQ. 

Start from Eq. 1)1(3 1) . and consider the following series 
expansions around x = 



erf (/ix) 



(Al) 



n=0 



2 (-1)" 
7^(2n-|- l)n!" 



(A2) 



This last series has an infinite radius of convergence for 
any finite /i. 

The complicated operator £^ can be also expanded in 
powers of x around a; = 0. Its expansion will only contain 
even powers of x because the hamiltonian of Eq. 1)15)) 
is even in x = vyi. As expected, the expansion of f 
does not play any role, so we do not consider its non- 
spherical components (moreover, in the end we are only 
interested in the spherically averaged pair density). The 
only important requirement is that remains bounded 
when /i — > 00 and x — *■ 0, as it happens for the Coulomb 
interaction 0,0]. We thus write 



(A3) 



fe=0 



By inserting Eqs. l)Aip . l)A2p and l)A3p into Eq. JTBl we 
find that the a„(/i) with odd n are zero (as expected 
from the fact that erf(/.ta;)/a; is an even function of x), 
while the even n coefficients with n > 2 diverge as /i 
increases, and, to leading order when \x — > 00, they are 
all proportional to ao{fi), 



a2fe+2(M) = ao(/i) 



,2k 



(2fc + 2)(2A: + 3) (2A: + 2)(2fc + 3) 



This relation shows that ijj'^{x) has the structure 
^P^^{x) = ao{^j.)[l + xpi{fix) + x'^p2{y.x) + ...] 



where 
Piiy) 



E 

00 



bkV 



2k+l 



P2{y) = 



(2fc + 2)(2A: + 3) 



,,2k 



^ I (2fc + 2)(2A: + 3) ^ 2i{2i + 1) 



(A4) 

(A5) 

(A6) 
I (AT) 



Equation l)A6p gives exactly the same function pi{y) of 
Eq. 1)22)) , and Eq. I)A4)I confirms that all the terms beyond 
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the ones considered in Eq. (|25|l contribute to Eq. H26|) to 
orders or higher. 

We can now clearly see where the discrepancy with the 
result of Ref. 10 comes from: the small-ri2 expansion of 
f^{fi2) only contains even powers of ri2 for any finite ^. 
It is thus incorrect to insert the odd coefficient ci of the 
Coulombic system in Eq. (|14|l . as it was done in Ref. fTol 
What happens, instead, is that all the even coefficients 
of the small-ri2 expansion of f^{r 12) diverge for large /i 
and they all contribute to the term oc in Eq. 
Furthermore, in Ref. |2 it was assumed that the leading 
order in the large-^ expansion of the on-top value /^(O) 
is 1/ IJ? , while we have shown that the correction to /^(O) 
with respect to the Coulombic case is of order 1///. 



APPENDIX B: THE CASE I = 1 

We insert the expansion of u^^{y) for large fi of Eq. I|18|l 
into the equivalent of Eq. 11711 for the case £ ^ 1, 



1 erf(y) 



M y 



= — £''u^(y). (Bl) 



The condition that the left-hand side be of order l/fi^ 
yields 



u^°°\y) = by' (B2) 
u(-i)(y) = byeviiy). (B3) 



By solving Eq. ljB3|) we find that the intracule f^{r 12) 
has, for large /x, the small-ri2 expansion 



r{ri2) 



/"(O) 



2 Bl 

1 + 2ri2qi{^iri2) + „ _ H 



where the function qi{y) is equal to 



91(2/) 



{2y' - 1) 



+ 



(B4) 



erf(y)(V + i;^ 
16^^^^) 



and Bl is a constant of integration that is not determined 
by the requirement that vanishes at ri2 = 0. The 
function qi{y) has the asymptotic behaviors 



91 (y - 



0) 



00 = 

' 4 30Fy 



(B6) 
(B7) 



Again, we see that if we fix ri2 = ^ 1, and then 
let — > 00 we find f^{ro) c>c rQ[l + ro/2 + ...], which 
is the parallel-spin cusp condition for the Coulomb in- 
teraction jl9j . But for any finite fi we have, for small 
7-12, F(ri2) oc rl2[l + rl22^l/5y^+ ■■■]■ The proof that 
Bl = 0, and thus of Eqs. H34|l and (|35|l is then completely 
analogous to the one for the case i — 0. 
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